The endpoint of the rp process on accreting neutron stars 
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We calculate the rapid proton (rp) capture process of hydrogen burning on the surface of an 
accreting neutron star with an updated reaction network that extends up to Xe, far beyond previous 
work. In both steady-state nuclear burning appropriate for rapidly accreting neutron stars (such as 
the magnetic polar caps of accreting X-ray pulsars) and unstable burning of Type I X-ray bursts, 
we find that the rp process ends in a closed SnSbTe cycle. This prevents the synthesis of elements 
heavier than Te and has important consequences for X-ray burst profiles, the composition of accreting 
neutron stars, and potentially galactic nucleosynthesis of light p nuclei. 
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I. INTRODUCTION 



Observations of X-ray binaries are an important source of information about neutron stars. In these systems a 
neutron star accretes hydrogen and helium rich matter from the envelope of a close companion star. Gravitational 
, energy is released when the matter reaches the neutron star and powers most of the observable X-ray radiation. 
' The accreted material is compressed, heated, and eventually undergoes thermonuclear burning. The released nuclear 
' energy, though small relative to the gravitational energy, is observable as Type I X-ray bursts at accretion rates 
. ^ 10~^Mo yr^*. For these rates, the nuclear burning proceeds unstably repeating on typical timescales of 

hours or days. In bright X-ray sources and in magnetically accreting X-ray pulsars the accretion rates reach or locally 
, exceed the Eddington limit. The accreted layer is then thermally stable and the thermonuclear burning proceeds in 
Ph' steady state ||J|] . 

Q When hydrogen is present, the nuclear burning proceeds via the rapid proton capture process (rp process) 0, a 

^ • sequence of proton captures and (3 decays responsible for the burning of hydrogen into heavier elements. We show 
' here, for the first time, that there is a natural termination of this process due to a closed cycle in the Sn-Te region 
C3 " (SnSbTe cycle). 

^ Understanding the rp process and especially its endpoint is a prerequisite for the interpretation of X-ray burst light 

• i-H . curves, which show a wide variety of characteristics including millisecond oscillations H] . Furthermore, the ashes of the 
nuclear burning replace the original crust of the neutron star and determine its thermal and electrical conductivities. 
The nuclear burning products thus have direct repercussions for the evolution of magnetic fields the quiescent 
luminosity in transient X-ray binaries [p^ (which could be used to discriminate against systems with accreting black 
holes), and the emission of potentially detectable gravitational waves from a rapidly rotating neutron star with a 
deformed crust 12|. There has also been an interest in the nuclear burning processes during X-ray bursts as a 



potential site for the nucleosynthesis of some light p nuclei like ^^^^^Mo and ^^'^^Ru that are severely underproduced 
in standard p process scenarios (for example |13|| ). This would require the escape of some of the burned material into 
the interstellar medium. 

Many early rp process simulations for X-ray bursts were based on reaction networks that ended at ^^Ni. For 
these calculations, it was argued that further proton captures were negligible for hydrogen consumption and energ 
production. However, several later studies were performed with larger networks ending in the Kr-Y region \ 
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or using a simplified 16 nuclei network up to Cd | [18[ . Recently a parameter study was carried out with an updated 
and extended network up to Sn All of these studies found that a significant fraction of the rp process reached 
the end of their respective networks. An rp process reaching the end of a reaction network up to Sn was also clearly 
demonstrated for the steady state burning regime for accretion rates between 20 and 40 times the Eddington limit 
(mEdd =8.8 104 g/cmVs) [|. 
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II. CALCULATIONS 



In this paper we explore the rp process beyond Sn in explosive and steady state burning to determine for the first 
time the endpoint of the rp process. Our one-zone X-ray burst model is based on the physics outlined in reference [pO| . 
Realistic ignition conditions have been calculated as a function of accretion rate and metallicity . For the X-ray 
burst model we assume an accretion rate of 0.1 mEdd, a flux out of the crust of 0.15 MeV per accreted nucleon, and 
a metallicity Z = 10~^. An accreted layer with solar metallicity would produce a similar burst at a higher accretion 
rate of 0.3 rriEdd- The steady state model is the same as described in reference p]. For all models we assume accretion 
of hydrogen and helium in solar proportions, and a surface gravity of 1.9x lO^^cm/s^. 

The nuclear reaction network includes all proton rich nuclei from hydrogen to xenon and was updated relative 
to the data described in ||l^. The theoretical Hauser-Feshbach reaction rates have been recalculated with the new 
Hauser-Feshbach code NON-SMOKER A more detailed discussion of the nuclear physics input and the X-ray 
burst model will be published in a forthcoming paper. 

III. RESULTS 

Figure |l| shows our results for the time integrated reaction flow during an X-ray burst. Ignition takes place at a 
density of 1.1 x 10^ gjcvc? and the burst reaches a peak temperature of 1.9 GK, with a rise time scale of «4 s, and a 
cooling phase lasting Ri200 s. Helium burns via the 3a reaction, and the ap process 0, a sequence of alternating (a,p) 
and (p,7) reactions into the Sc region. These helium burning processes provide the seed nuclei for the rp process. 
The rp process reaction flow reaches the Sn isotopes in the ^^Sn-^^^Sn range « 80 s (time for half maximum) after 
the burst peak and proceeds then along the Sn isotopic chain towards more stable nuclei. 

Processing beyond Sn occurs if the corresponding Sb isotone is sufficiently proton bound for the (7,p) photodisinte- 
gration to be small. This occurs at ^"^Sn. However, after two proton captures a strong ^°^Te(7,a) photodisintegration 
rate cycles the reaction flow back to ^''■^Sn. The reaction path is characterized by a cyclic flow pattern, the SnSbTe 
cycle which represents the endpoint for the rp-process reaction flow towards higher masses (see figure 0). The SnSbTe 
cycle forms because the neutron deficient i06-i08rpg isotopes are a unbound by «4 MeV. In fact, ^^Te is a known 
ground state a emitter j2^]. A fraction of the reaction flow proceeds via /3 decay of ^°^Sn into ^°^Sn, and the reaction 
sequence ^°^Sn(p,7)-^°''Sb(p,7)^''^Te(7,a)^°'*Sn leads to a second, weaker cycle. Calculations with different ignition 
conditions confirm that the rp process cannot proceed beyond the SnSbTe cycles. 

A previous calculation of the rp process in steady state burning found that most material accumulated at the end of 
the network (the Sn isotopes) for an accretion rate of 40mEdd [|| ■ Figure |l| shows the reaction flow at that accretion 
rate. We find that the rp process ends in a similar SnSbTe cycle as in X-ray bursts. Some of the material is now cycled 
back via ^''^Sb(p,a), which successfully competes with ^°^Sb(p,7) at steady state burning conditions. Calculations 
at different accretion rates show that the rp process can never overcome the closed SnSbTe cycle. For steady state 
burning, we are now able to compute accurately the composition of the ashes for all accretion rates. 

The SnSbTe cycle impacts the light-curve of X-ray bursts and the consumption of hydrogen. This is illustrated 
in figure H, which shows the correlation between the X-ray burst luminosity, the abundances of some important long 
lived nuclei (waiting points) in the rp process, and the hydrogen and helium abundances. Clearly the slow hydrogen 
burning via the rp process beyond ^^Ni is responsible for the extended burst tail. The SnSbTe cycle builds up the 
abundance of the longest lived nucleus in the cycle, ^°'*Sn (20.8 s half-life), and produces helium towards the end of 
the burst. This triggers an increase in the 3a flow and subsequently an increase in energy production and hydrogen 
consumption. As a consequence, the burst lasts longer and hydrogen is completely burned. 

The SnSbTe cycle also affects the composition of the rp process ashes, shown for the X-ray burst and the steady 
state calculation in figure ^. The limitation imposed on the rp process by the SnSbTe cycle is clearly reflected in the 
lack of nuclei heavier than A — 107. Nevertheless we obtain a broad distribution of nuclei in the A — 64-107 mass 
range. This is a result of the long lived waiting point nuclei along the rp process reaction path which store some 
material until the burning is over. The late helium production in the SnSbTe cycle broadens this distribution further. 

IV. CONCLUSIONS 

To summarize, we have shown that the synthesis of heavy nuclei via the rp process is limited to nuclei with 
Z < 54 due to our newly discovered SnSbTe cycle. The existence of a SnSbTe cycle under all rp process conditions 
is a consequence of the low, experimentally known |Q a separation energies of the I06,i07,i08,i09rpg isotopes and is 
therefore not subject to nuclear physics uncertainties. However, because of the uncertainties in the proton separation 
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energies of the Sb isotopes there is some uncertainty in the relative strength of the SnSbTe sub-cycles closed by (7, a) 
photodisintegration on ^°^Te, ^°^Te, and ^°®Te. This will be discussed in a forthcoming paper. 

A likely consequence of the SnSbTe cycle for accreting neutron stars is that the matter entering the crust is composed 
of nuclei lighter than A «107. The only way to bypass the SnSbTe cycle would be a pulsed rp process, where between 
pulses matter could decay back to stable nuclei. This could happen during so called dwarf bursts, which have been 
suggested to be secondary bursts produced by reignition of the ashes |^5| . However, this would require some unburned 
hydrogen in the burst ashes (see discussion below) or extensive vertical mixing p4[ |. 

Our calculations give a strong indication that the synthesis of nuclei beyond °^i and especially into the ^ = 100 
mass region in hydrogen rich bursts leads to extended energy production. This might explain the long duration (100 
seconds) bursts seen from, for example, GS 1826-24 |p6|| . 

The importance of the endpoint of the rp process for the consumption of hydrogen in X-ray bursts has been discussed 
extensively. Previous calculations limited to nuclei up to ^^Ni always found unburned hydrogen. Later burning of 
this hydrogen was proposed to explain short burst intervals |25| , occasionally observed extended X-ray flares , and 
spin down in miUisecond oscillations, as recently discovered during the tail of a long X-ray burst from 4U 1636-53 [ p8| . 
It was pointed out before that depending on the endpoint of the rp process these models might largely overestimate 
the amount of unburned hydrogen [T^ . Our calculations indicate indeed that the SnSbTe cycle reduces the amount 
of unburned hydrogen substantially. Within the framework of the one-zone model we find that all the hydrogen is 
burned, a full multizone hydrodynamic burst calculation is needed to confirm this result. 

In our X-ray burst calculation, we find large overproduction factors (produced abundance to solar abundance « 10^) 
for the p nuclei ^^Ru, ^°^Pd, and ^°^Cd. Overproduction factors of this order have been suggested to be sufficient to 
explain the origin of the solar system abundances of these nuclei if « 1% of the burned material is somehow ejected 
p9| . However, in order to accurately calculate the composition of the ejected material the ejection mechanism has to 
be identified. Furthermore we find that the synthesis of p nuclei is most likely associated with X-ray bursts having 
long tails. As these bursts are rare, we conclude that only a small fraction of X-ray bursts are likely to produce light 
p nuclei in the A ~ 90 — 106 mass range. Therefore, considerably larger amounts of matter would have to be ejected 
in X-ray bursts to account for the galactic abundances of light p nuclei than previously assumed. Nevertheless it is 
striking that the SnSbTe cycle happens to naturally limit rp process nucleosynthesis to light p nuclei. 

The authors would like to thank R. Boyd and E. Roeckl for helpful discussions. Part of this work was carried 
out under NSF contract PHY-95-28844, PHY-99-01133 and PHY-99-07949, under NASA grant NAG5-8658, and 
under and the Swiss NSF contracts 2124-055832.98 and 2000-061822.00. L.B. is a CottreU Scholar of the Research 
Corporation. V.B. was supported through the Graduate School of the University of Notre Dame. 
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FIG. 1. The time integrated reaction flow above Ga during an X-ray burst and for steady state burning. Shown are reaction 
flows of more than 10% (solid line) and of 1-10% (dashed line) of the reaction flow through the 3a reaction. 
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FIG. 2. The reactions in the SnSbTe cycles during an X-ray burst. In the case of proton captures the arrows indicate the 
direction of the net flow, the difference of the flow via proton capture and the reverse flow via (7,p) photodisintegration. The 
line styles are the same as in figure ^. 
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FIG. 3. Luminosity, nuclear energy generation rate, and the abundances of hydrogen, helium and the important waiting 
point nuclei as functions of time during an X-ray burst. For comparison, the nuclear energy generation rate is also shown as a 
dashed line together with the luminosity, though it is out of scale during the peak of the burst. The mass of the accreted layer 
is 4.9x10^^ g. 
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FIG. 4. Final abundance distribution as functions of mass number for an X-ray burst, and for steady state burning at an 
accretion rate of 40mEdd- 
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